Properties of Elastic Waves in a non-Newtonian (Maxwell) 
Fluid-Saturated Porous Medium 

^^ ' David Tsiklauri (tsikd@astro.warwick.ac.uk) 

r\] , Physics Department, University of Warwick, Coventry, CV4 7AL, UK 

j^ \ Igor Beresnev (beresnevOiastate . edu) 

Department of Geological and Atmospheric Sciences, Iowa State University, 253 
Science I, Ames, IA 50011-3212, U.S.A. 

00 ' 

CN) ■ February 2, 2008 

Abstract. The present study investigates novelties brought about into the classic 
Biot's theory of propagation of elastic waves in a fluid-saturated porous solid by 

0-, inclusion of non-Newtonian effects that are important, for example, for hydrocar- 

bons. Based on our previous results (Tsiklauri & Beresnev: 200f, Phys. Rev. E, 63, 
046304), we have investigated the propagation of rotational and dilatational elastic 

j^JT)' waves, through calculating their phase velocities and attenuation coefficients as a 

function of frequency. We found that the replacement of an ordinary Newtonian 
fluid by a Maxwell fluid in the fluid-saturated porous solid results in: (a) an overall 
increase of the phase velocities of both the rotational and dilatational waves. With 
the increase of frequency these quantities tend to a fixed, higher, as compared to 
the Newtonian limiting case, level which is not changing with the decrease of the 
Deborah number a. (b) the overall decrease of the attenuation coefficients of both 
the rotational and dilatational waves. With the increase of frequency these quantities 
tend to a progressively lower, as compared to the Newtonian limiting case, levels as 
a decreases, (c) Appearance of oscillations in all physical quantities in the deeply 
non-Newtonian regime. 
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c/3 ' 1. Introduction 

O 

• i— i 

Apart from fundamental interest, there are at least three major rea- 
sons to study the dynamics of fluid in porous media under oscillatory 
Q_c pressure gradient and oscillating pore walls, as well as to investigate 

propagation of elastic waves in porous media. 

First, in petroleum geophysics, regional exploration seismology needs 
direct methods of discovering oil-filled bodies of rock, and these should 
be based on models of propagation of elastic waves in porous media 
with realistic fluid Theologies (Carcione h. Quiroga-Goode, 1996). 

Second, the investigation of the dynamics of fluid in porous me- 
dia under oscillatory pressure gradient is of prime importance for the 
recently emerged technology of acoustic stimulation of oil reservoirs 
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(Beresnev & Johnson, 1994), (Drake & Beresnev, 1999). It is known 
that the natural pressure in an oil reservoir generally yields no more 
than approximately 10 percent oil recovery. The residual oil is difficult 
to produce due to its naturally low mobility, and the enhanced oil 
recovery operations are used to increase production. It has been exper- 
imentally proven that there is a substantial increase in the net fluid 
flow through porous space if the latter is treated with elastic waves 
(Beresnev & Johnson, 1994), (Drake & Beresnev, 1999). 

Third, in the environment conservation, treatment of ground water 
aquifers contaminated by organic liquids, such as hydrocarbons, by 
elastic waves proved to be successful for quick and efficient clean up 
(Beresnev & Johnson, 1994), (Drake & Beresnev, 1999). 

A quantitative theory of propagation of elastic waves in a fluid- 
saturated porous solid was formulated in the classic papers by Biot 
(Biot, 1956a, b). One of the major findings of Biot's work was that 
there is a breakdown in Poisseuille flow above a certain characteristic 
frequency specific to the fluid-saturated porous material. Biot theoret- 
ically studied this phenomenon by considering the flow of a viscous 
fluid in a tube with longitudinally oscillating walls under an oscillatory 
pressure gradient. In Biot's theory, the two-phase material is considered 
as a continuum and the microscopic, pore-level effects are ignored. As a 
reminder, the theory assumes that: (a) the wavelength is large with re- 
spect to the dimensions of pores in order to make continuum mechanics 
applicable; this also implies that scaterring dissipation is negligible; (b) 
the displacements are small, therefore the macroscopic strain tensor is 
related to them by the lowest second-order approximation; (c) the liquid 
phase is continuous, such that the pores are connected and isolated 
pores are treated as part of solid matrix; and (d) the permeability is 
isotropic and the medium is fully saturated. 

Biot demonstrated the existence of the two kinds of compressional 
waves in a fluid-saturated porous medium: the fast wave for which the 
solid and fluid displacements are in phase, and the slow wave for which 
the displacements are out of phase. At low frequencies, the medium 
does not support the slow wave as it becomes diffusive. On the other 
hand, at high frequencies, tangential slip takes place, inertial effects 
dominate, and the slow wave becomes activated. 

Biot's theory can be used to describe interaction of fluid-saturated 
solid with the sound for a classic Newtonian fluid; however, oil and 
other hydrocarbons often exhibit significant non-Newtonian behavior. 
In paper I (Tsiklauri & Beresnev, 2001), we have incorporated non- 
Newtonian effects into the classical theory of Biot (Biot, 1956a, b). 
Using the recent results of del Rio et al. (1998), who presented a study of 
enhancement in the dynamic response of a viscoelastic (Maxwell) fluid 
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flowing through a stationary (non-oscillating) tube under the effect of 
an oscillatory pressure gradient, we have combined their theory with 
the effect of the acoustic oscillations of the walls of the tube introduced 
by Biot (Biot, 1956a,b), thus providing a complete description of the 
interaction of Maxwell fluid, filling the pores, with acoustic waves. We 
have generalized the expression for the function F(k), which measures 
the deviation from Poisseuille flow friction as a function of frequency 
parameter k (Tsiklauri & Beresnev, 2001). As a next step, in the 
present work we investigate the propagation of rotational and dilata- 
tional elastic waves through the porous medium filled with Maxwell 
fluid, by calculating their phase velocities and attenuation coefficients 
as a function of frequency. 

This paper is organized as follows: we formulate theoretical basis for 
our numerical calculations in section 2. In sections 3 and 4 we study 
numerically properties of the rotational and dilatational elastic waves, 
respectively, and, finally, in Section 5 we close with a discussion of our 
main results. 



2. Theory 

The theory of propagation of elastic waves in a fluid-saturated porous 
solid was formulated by Biot (Biot, 1956a, b). He demonstrated that the 
general equations which govern propagation of rotational and dilata- 
tional high-frequency waves in a fluid-saturated porous medium are the 
same as in the low-frequency range provided the viscosity is replaced 
by its effective value as a function of frequency. In practice, it means 
replacing the resistance coefficient b by bF(n). 

The equations describing dynamics of the rotational waves are (Biot, 
1956a,b) 

3j2 (Pud + Pnfi) + bF ^g- t ^ ~ ") = NV ^ C 1 ) 

-Q^{pi2^ + P22^)-bF{K)- (£3-fi) = 0, (2) 

where, pu,p\2 and P22 are mass density parameters for the solid and 
fluid and their inertia coupling; uj = curl u and Q = curl U describe 
rotations of solid and fluid with u and U being their displacement 
vectors, while the rigidity of the solid is represented by the modulus 
N. Substitution of a plane rotational wave of the form 



u 



Cie i(fa+xt) ) n = C 2 e i( - lx+xt) , (3) 
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into Eqs.(l) and (2) allows us to obtain a characteristic equation 

Nl 2 

— = E r - iEi, (4) 

PX 2 

where I is wavenumber, x = 2nf is wave cyclic frequency, p = p\\ + 
2pi2 + P22 is the mass density of the bulk material. 

The real and imaginary parts of Eq. (4) can be written as 



(711722 - 71 2) (722 + €2) +722 62 + e? + ej 
(722 + e 2 ) 2 + 4 



7-1 v /11 rzz IV2JK 1^^ ' ^^J ' l^^-t ' "-1 1 c 2 /r\ 

^ — 77 , , \9 , ,2 ' W 

and 



^ = t 1 7^ ( 6 ) 

(722 + e 2 ) 2 + ef 



£1(712 + 722 ^ 2 

(722 + £2) 2 + 

where 7^ = p^/p, e x = (71 2 + 7 22) (fc/f) Re[F(n)] = 
(7i2 + 722)(fc/f)Re[F(5y/f/f c )}, e 2 = (712 + 722) (f c / f)Im[F(n)\ = 
(712 +722) (fc/f) I rn [F(5^f/f c )]. The function F(k) was written here 
more conveniently as a function of frequency / , i.e. F(k) = F(5y/f/f c ) 
(Biot, 1956a, b), where 5 is a factor dependent on pore geometry. For 
the hollow cylinder- like pores, 5 = y/8 (Biot, 1956a, b) and we use this 
value throughout the paper. f c is the critical frequency above which the 
Poisseuille flow breaks down, and it equals b/(2np2) = b/(2irp(ji2 + 

722))- 

In order to obtain phase velocity and attenuation coefficient of the 
rotational waves, we put I = Re[l] + ilm[l]. Thus, the phase velocity 
is then v r = x/|ite[Z]|. Introducing a reference velocity as V r = y/N/p, 
we obtain the dimensionless phase velocity as 

v jl = yl (7) 



E 2 + E 2 + E, 



r 



1/2- 



To obtain the attenuation coefficient of the rotational waves, we 
introduce a reference length, L r , defined as L r = V r /(2irf c ). The length 
x a represents the distance over which the rotational wave amplitude is 
attenuated by a factor of 1/e. Therefore we can construct the dimen- 
sionless attenuation coefficient as L r /x a , 

1/2 



L r f 



E 2 + E 2 - E r 



The equations describing dynamics of the dilatational waves are 
(Biot, 1956a,b) 

V 2 (Pe + Qe) = QpiPiie + Pi 2 e) + bF ( K )^( e " e )> ( 9 ) 



tipml765.tex; 2/02/2008; 1:52; p. 4 



Elastic Waves in a Maxwell fluid-Saturated matrix 5 

V 2 (Qe + Re) = Q^(pi2e + p 22 e) - *>F( K )-^(e " e )> ( 10 ) 

where, P, Q and R are the elastic coefficients, e = div u and e = div J7 
are the divergence of solid and fluid displacements. Again, substitution 
of a plane dilatational wave of the form 

e = de i{lx+xt \ e = C 2 e i{lx+Xt \ (11) 

into Eqs.(9) and (10) allows us to obtain a characteristic equation 

(z- Zl ){z-z 2 ) + iM(z-l) = 0, (12) 

where z = l 2 V 2 /x 2 , V 2 = (P + R + 2Q)/p represents the velocity of 
a dilatational wave when the relative motion between fluid and solid 
is absent, z± :2 = V 2 /V 2 2 with V± :2 being the velocities of the purely- 
elastic waves with subscripts 1,2 referring to the two roots of Eq.(12), 
and finally M = (ei + ie 2 )/(a u cr 22 - cr 2 2 ) with an = P/{P + R + 2Q), 
a 22 = R/(P + R + 2Q) and a 12 = Q/(P + R + 2Q). 

Eq.(12) has two complex roots zj and zn- Phase velocities of the 
two kinds of dilatational waves can be defined as 

1 '"" l (13) 



V c Re[^zJY v c Re[^zny 

while the corresponding attenuation coefficients can be also introduced 
as 

h. = Im [ y /^\f i^ = j m [ v ^ 7 ]Z. ( M ) 

XI Jc XII Jc 

In paper I, we generalized Biot's expression for F(k) to the case of 
a non-Newtonian (Maxwell) fluid, which reads 



i K^Ji + n 2 /a J\{nyi + k 2 / 'a) / 'Jo(ny/ 'i + k 2 /a) 
F(k)-- 



4 (1 - ik 2 /a) 

-l 

(15) 



2Ji(/c v /« + K 2 /a) 



K\Ji + k 2 J 'aJo{n-\/ i + K 2 /a) 



Here, k = a^Juj/u is the frequency parameter, a is the radius of the 
pore, v = 7]/p is the ratio of the viscosity coefficient to the fluid mass 
density, Jo and J\ are the Bessel functions, and, finally, a denotes the 
Deborah number del Rio et al. (1998), which is defined as the ratio of 
the characteristic time of viscous effects t v = a 2 jv to the relaxation 
time t m , i.e., a = t v /t m = a 2 /(vt m ). 

Eq.(15) was derived by solving the equations of incompressible hy- 
drodynamics, namely, the continuity equation, linearized momentum 
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equation, and rheological equation of a Maxwell fluid, in the frequency 
domain for a cylindrical tube whose walls oscillate harmonically in time. 
By calculating the ratio of the total friction force exerted on the tube 
wall to the average velocity of a Maxwell fluid, and noting that F{n) is 
proportional to this ratio, we generalized the classical result obtained 
by Biot (see details in Paper I). 

As noted in del Rio et al. (1998), the value of the parameter a 
determines in which regime the system resides. Beyond a certain crit- 
ical value (a c = 11.64), the system is dissipative, and viscous effects 
dominate. On the other hand, for small a (a < a c ), the system exhibits 
viscoelastic behavior which we call the non-Newtonian regime. Note, 
that the Newtonian flow regime can be easily recovered from Eq.(15) 
by putting a — > oo. 

In order to investigate the novelties brought about into classical 
Biot's theory of propagation of elastic waves in porous medium (Biot, 
1956a,b) by the inclusion of non-Newtonian effects, we have studied the 
full parameter space of the problem. We have calculated the normalized 
phase velocities and attenuation coefficients for both rotational and 
dilatational waves using our general expression for F(k) given by Eq. 
(15). 



3. Numerical Results for Propagation of Rotational Waves 

In all our numerical calculations, we have used polynomial expansions 
of Jo and J\ with absolute error not exceeding 10~ 6 %. Thus, our 
calculation results are accurate to this order. Also, in order to catch a 
true oscillatory structure of our solutions (see below), number of data 
points in all our plots is 10000, as opposed to paper I where only 100 
points per curve were taken. 

In all forthcoming results, we calculate phase velocities and at- 
tenuation coefficients for the case 1 from Table I taken from (Biot, 
1956b), which is a n = 0.610, a 22 = 0.305, a 12 = 0.043, 711 = 0.500, 
722 = 0.500, 712 = 0, zi = 0.812, and z 2 = 1.674. 

We calculated normalized phase velocity of the plane rotational 
waves, v r /V r , and the attenuation coefficient L r /x a using our more 
general expression for F(n) (Maxwell fluid filling the pores) given by 
Eq. (15). 

In Fig. 1 we plot phase velocity v r /V r as a function of frequency 
for the three cases: the thick curve corresponds toa^oo (Newtonian 
limit), the dashed curve corresponds to a slightly sub-critical value of 
a = 10 (recall that a c = 11.64) and thin solid curve corresponds to 
the case of a = 1. Note that the a — ► 00 case perfectly reproduces the 
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curve 1 in Fig. 5 from (Biot, 1956b). For a = 10 we notice a deviation 
from the classic Newtonian behavior in the form of overall increase 
of phase velocity and appearance of small oscillations on the curve, 
which means that we have entered the non-Newtonian regime. Note 
that when a = 1 the phase velocity settles at somewhat higher value 
and this onset happens already for smaller frequencies than in the case 
of Newtonian fluid. Also, much more pronounced oscillatory structure 
of the solution can be observed. 

Fig. 2 shows the attenuation coefficient L r /x a of the rotational wave 
as a function of frequency for the three values of a: the thick curve cor- 
responds to a — > oo (Newtonian limit), the dashed curve corresponds to 
a slightly sub-critical value of a = 10 and thin solid curve corresponds 
to the case of a = 1. Note that a — ► oo case coincides with curve 1 in 
Fig. 6 from (Biot, 1956b). For a = 10, there is a noticeable deviation 
from the classic Newtonian behavior in the form of overall decrease of 
the attenuation coefficient and appearance of small oscillations on the 
curve indicating that the wave has entered the non-Newtonian regime. 
For a = 1, the attenuation coefficient settles at a somewhat lower value, 
and this happens already for smaller frequencies than in the case of 
Newtonian fluid. Also, much more pronounced oscillatory structure of 
the solution can be noticed. 



4. Numerical Results for Propagation of Dilatational Waves 

We calculated normalized phase velocities of the plane dilatational 
waves, vj/V c and vu/V c , and the attenuation coefficients L c /xi and 
L c /xn using our more general expression for F(k) (Maxwell fluid filling 
the pores) given by Eq. (15). 

In Fig. 3 we plot phase velocity vi/V c as a function of frequency for 
the case of a — > oo, in order to recover the Newtonian limit obtained 
by Biot. Note that this case reproduces curve 1 in Fig. 11 from (Biot, 
1956b). 

Fig. 4 shows phase velocity vi/V c as a function of frequency for 
the case of a = 1, corresponding to the deeply non-Newtonian regime. 
We notice again appearance of an oscillatory structure of the solution. 
Also, phase velocity vi/V c settles at a somewhat higher value than in 
the Newtonian case, and this happens already for smaller frequencies. 

In Fig. 5 we plot phase velocity vu/V c as a function of frequency 
for the cases of a — > oo and a = 1. Note that the case of a — > oo (thick 
solid line) perfectly reproduces the Newtonian limit obtained by Biot 
curve 1 in Fig. 12 from (Biot, 1956b). For a = 1, again we notice an 
oscillatory structure of the solution. Besides, we observe that the phase 
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velocity vjj/V c settles again at a somewhat higher value than in the 
Newtonian case. 

In Fig. 6 we plot the attenuation coefficient L c /xj as a function 
of frequency for the three cases: the thick curve corresponds to a — > 
oo (Newtonian limit), the dashed curve corresponds to a slightly sub- 
critical value of a = 10 and thin solid curve corresponds to the case 
of a = 1. The a — ► oo case reproduces curve 1 in Fig. 13 from (Biot, 
1956b). For a = 10, we notice deviation from the classic Newtonian 
behavior in the form of the overall increase of the attenuation coefficient 
and appearance of small oscillations on the curve, which indicates that 
the wave has entered the non-Newtonian regime. The large spike at 
low frequencies is also due to non-Newtonian effects. For the case of 
a = 1 the attenuation coefficient settles at somewhat lower values, 
and this happens already for smaller frequencies than in the case of a 
Newtonian fluid. Also, a much more pronounced oscillatory structure 
of the solution can be noticed. 

In Fig. 7 we plot the attenuation coefficient L c /xn as a function of 
frequency for the three cases: the thick curve corresponds to a — ► oo 
(Newtonian limit), the dashed curve corresponds to a = 10 and thin 
solid curve corresponds to the case of a = 1. Note that the a — > oo case 
perfectly matches curve 1 in Fig. 14 from (Biot, 1956b). For a = 10, 
we notice a deviation from the classic Newtonian behavior in the form 
of the overall decrease in the attenuation coefficient and appearance of 
small oscillations on the curve. The jump at f/f c = 1 (dashed curve) 
should be attributed to the non-Newtonian effects. For the case of a = 1 
the attenuation coefficient settles again at somewhat lower value and 
this happens already for smaller frequencies than in the case of the 
Newtonian fluid. Also, much more pronounced oscillatory structure of 
the solution can be noticed again. 



5. discussion 

In this paper, we have studied the non-Newtonian effects in the propa- 
gation of elastic waves in porous media by calculating phase velocities 
and attenuation coefficients of the rotational and dilatational waves 
as a function of frequency. Originally, Biot (Biot, 1956a, b) performed 
similar analysis for a Newtonian fluid-saturated porous medium. Us- 
ing our recent results [Paper I], and motivated by a current need in 
models of propagation of elastic waves in porous media with realistic 
fluid rheologies, we have generalized the work of Biot to the case of a 
non-Newtonian (Maxwell) fluid-saturated porous medium. 
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In summary, we found that replacement of an ordinary Newtonian 
fluid by a Maxwell fluid in the fluid-saturated porous medium results 
in 

— an overall increase of the phase velocities of both the rotational and 
dilatational waves. With the increase of frequency these quantities 
tend to a fixed, higher, as compared to the Newtonian limiting 
case, level which is not changing with the decrease of the Deborah 
number a. 

— the overall decrease of the attenuation coefficients of both the 
rotational and dilatational waves. With the increase of frequency 
these quantities tend to a progressively lower, as compared to the 
Newtonian limiting case, levels as a decreases. 

— Appearance of oscillations in all physical quantities in the deeply 
non-Newtonian regime when a <C a c = 11.64. 

The investigation of properties of elastic waves is important for a 
number of applications. The knowledge of phase velocities and attenu- 
ation coefficients of elastic waves in a realistic [such as saturated with 
Maxwell fluid] porous medium is necessary, for example, to guide the 
oil- field exploration applications, acoustic stimulation of oil producing 
fields (in order to increase the amount of recovered residual oil), and 
the acoustic clean up of contaminated aquifers (Beresnev & Johnson, 
1994), (Drake & Beresnev, 1999). 

The idea of the paper was to use the function, F(k), that measures 
the deviation from Poisseuille flow friction, extended to Maxwell flu- 
ids, and to substitute it into Biot's equations of poroelasticity without 
changing the latter. However, Biot's equations have been derived under 
a number of assumptions. One of these assumptions is that deviatoric 
(shear) components of the macroscopic stress in the fluid are negligible 
Pride et al. (1992). Pride et al. (1992) have shown that this assumption 
is justified when ijuj <C G, where r\ is the dynamic viscosity of the 
fluid, to is frequency, and G is frame shear modulus. Simple analysis 
shows that for typical Newtonian fluids such as water, this condition 
is only violated at frequencies u> > 10 9 1/s, or / = uj/(2tt) > 10 8 Hz. 
Thus, for all frequencies below 1 MHz Biot's assumption is justified. 
However, when we introduce the Maxwell fluid, the situation changes 
in that we introduce the real (in addition to imaginary) shear stresses. 
In summary, for any rheology (including Maxwellian) Biot's theory is 
only valid if macroscopic shear stresses are negligible. In order to prove 
that, we note from the rheological equation for a Maxwell fluid 

df „^ _ 

tm-^7 = -rjS/V - T, 
Ot 
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where f represents the viscous stress tensor, that in the frequency 
domain we can effectively obtain 

f = —rf\7v/{l + it m Lu). 

This means that we can roughly replace rj in all our estimates with 
r//(l + it m (jj). There are two limiting cases. When uj <C l/t m , then the 
fluid is effectively Newtonian and Biot's theory is valid. When to S> 
l/t m , i.e. when the fluid is essentially non-Newtonian, we effectively 
have rj' = rj/(it m uj), which in this case is smaller than 77 in absolute 
value. Thus, when substituted into the shear stress, S, it produces 
S = ir/'u = rj/t m , which is smaller than r/uj. Therefore, we conclude 
that inequality rj'uo <C G still holds for the Maxwellian fluid, i.e. Biot's 
equations are valid for Maxwell rheology. 

This study, similarly to the results of Paper I, has clearly shown the 
transition from dissipative to non-Newtonian regime in which sharp 
oscillations of all physical quantities are found. We would like to com- 
ment on these unexpected strong oscillations that were demonstrated 
by our numerical analysis. The results are based on equation (15), 
which has been derived for a circular cylindrical geometry. This is the 
same geometry that was used in classical works of Biot and others. 
For Newtonian fluids the use of such an idealized geometry for porous 
materials was backed up by an analysis that showed that the results 
are not very sensitive to the particular geometry (see, e.g., Johnson et 
al. (1987)). Of course, the magnitude of these oscillations depends on 
fluid parameters and permeability, and may not be as high for many 
fluids. However, even if parameters are such that oscillations are large, 
it is unclear at this stage whether this oscillatory behavior will hold 
for more realistic geometry, i.e. when curved pore walls (tortuosity) are 
considered. There is a possibility that with tortuosity effects included 
the obtained oscillations will be smeared. However, our goal was to 
constrain ourselves with simple cylindrical geometry, and a separate 
study is needed to analyze the tortuosity effects on our results. 
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Figure captions 

Fig. 1 Behavior of dimensionless, normalized phase velocity of the 
rotational wave, v r /V r , as a function of frequency. The thick solid curve 
corresponds to the Newtonian limit when a — ► oo, while dashed and 
thin solid curves represent the non-Newtonian cases a = 10 and a = 1 
respectively. 

Fig. 2 Behavior of dimensionless, normalized attenuation coefficient 
of the rotational wave, L r /x a , as a function of frequency. The thick solid 
curve corresponds to the Newtonian limit when a — ► oo, while dashed 
and thin solid curves represent the non-Newtonian cases a = 10 and 
a = 1 respectively. 

Fig. 3 Behavior of dimensionless, normalized phase velocity of the 
dilatational wave, vi/V c , as a function of frequency. Here the Newtonian 
limiting case, when a — ► oo, is plotted. 

Fig. 4 Same as in Fig. 3 but for a = 1. 

Fig. 5 Behavior of dimensionless, normalized phase velocity of the 
dilatational wave, vu/V c , as a function of frequency. Here the Newto- 
nian limiting case, when a — > oo, is plotted with thick curve, while the 
thin one corresponds to the non-Newtonian case of a = 1. 

Fig. 6 Behavior of dimensionless, normalized attenuation coefficient 
of the dilatational wave, L c /xj, as a function of frequency. The thick 
solid curve corresponds to the Newtonian limit when a — ► oo, while 
dashed and thin solid curves represent the non-Newtonian cases a = 10 
and a = 1 respectively. 

Fig. 7 Behavior of dimensionless, normalized attenuation coefficient 
of the dilatational wave, L c /xu, as a function of frequency. The thick 
solid curve corresponds to the Newtonian limit when a — > oo, while 
dashed and thin solid curves represent the non-Newtonian cases a = 10 
and a = 1 respectively. 
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